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Abstract. 

We present new observations of the famous starburst galaxy NGC 5253 which owes its celebrity to possibly being the 
youngest and closest starburst galaxy known. Our observations in the infrared and millimeter contribute to shed light on the 
properties of this interesting object. We have used our new data along with data from the literature to study the properties of the 
young stellar clusters present in NGC 5253. We find that the brightest optical clusters are all characterized by a near-infrared 
excess that is explained by the combined effect of extinction and emission by dust. For the brightest infrared cluster we model 
the spectral energy distribution from the optical to the radio. We find that this cluster dominates the galaxy emission longward 
of 3 yum, that it has a bolometric luminosity of 1.2 10' L Q and a mass of 1.2 10 6 M Q , giving L/M ss 10 3 . The cluster is obscured 
by 7 mag of optical extinction produced by about 1.5 10 5 M Q of dust. The dust properties are peculiar with respect to the dust 
properties in the solar neighbourhood with a composition characterized by a lack of silicates and a flatter size distribution 
than the standard one, i.e. a bias toward larger grains. We find that NGC 5253 is a striking example of a galaxy where the 
infrared-submillimeter and ultraviolet-optical emissions originate in totally decoupled regions of vastly different physical sizes. 
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£j ; 1. Introduction 

^ It has been found that very young clusters can be deeply embedded in dusty cocoons even when hosted in extremely low metal- 
kJ ' licity environments. The most prominent example of this class o f objects is certainly the Super-Star Cluster (SSC) observed in the 
• i-h ! extremely metal deficient blue compact galaxy SBS 0335-052 jThuan et alJl999tlHunt et all200ll) . Modeling the SED of this ob- 
' ject indicates that the radiation from a few /mi to the far-infrared (FIR) comes from a luminous (> 10 9 L ), massive (> lO 6 !^) 
^ \ star c luster embedded in a dust cocoon where large grains dominate, contrary to interstellar medium dust jPlante & Sauvagd 
- - ■ I2002I) . The amount of dust derived by fitting the SED, about 10 5 M , is surprisingly high for a very low metallicity environ- 
ment and it opens several interesting questions about the star formation history of the galaxy and the processes involved in dust 
formation. Clusters with similar properties and possi bly even m ore embedded, since they are not detected at ~2 /vm, have been 
observed in the blue compact galaxy He 2-10 JVacca et alJl2002h . It is found that these clusters have properties very similar to 
those of galactic ultra dense HII regions and that they are probably the main source of the observed far-infrared flux. In SBS 0335- 
052, the dust-enshrouded SSC contains about 10 times more 07V stars than its visible counterparts (Plante & Sauvage 2002). 
In He 2-10, the emb edded source represent > l/9 th of the total UV-emitting star content of the galaxy ( iKobulnickv & Johnson! 
ll999HConti & Vaccalll994. In the se two galaxies, the visib le super-star clusters account for almost all of the recent visible star 
formation activity dConti & Vaccal ll994: Thuan et a 1. 1997). Therefore these two examples show that a substantial, if not domi- 
nant, fraction of the current star formation can be made invisible in the optical-UV by dust obscuration even in low-metallicity 
environments. Were this to be a general property, it would obviously have important consequences on the interpretation of the 
high redshift universe where most star formation should occur in low-metallicity environments. 

In this paper, we study the case of NGC 5253, with the purpose of demonstrating that it is also a member of the class of 
galaxies that host dust-enshrouded SSCs and to strengthen the relevance of the phenomenon in the study of starburst objects. 

Send offprint requests to: L. Vanzi 

* Based on observations obtained at the ESO telescopes of La Silla and Paranal, program 69.B-0345; and on observations with ISO, an 
ESA project with instruments funded by ESA Member States (especially the PI countries: France, Germany, the Netherlands and the United 
Kingdom) and with the participation of ISAS and NASA. 
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The starburst galaxy NGC 5253 is a blue dwarf irregular galaxy which has gained increasing interest during the years because 
of a few important facts. (1) Among the starburst galaxies it is one of the closest e xamples. Its exac t distance is still a matter of 
debate, with the most recent values being either 3.3+0.3 Mpc jGibson etaPEo OO) or 4.0±0.3Mpc jThim et alJl2003l) . Here we 
arbitrarily choose to use a distance of 3.3 Mpc, for which 1 "is 16pc. (2) It is most likely the younge st starburst known with a large 
numb er of clusters detected whose ages are in the 2-50 Myr range (see e.g. Ivan den Bergh1ll980t iGoriian lfl"996t ICalzetti et alJ 
119971) . ( 3) Although its metal content is not extremely low, as often observed in dwarf galaxies, its metallicity is sub-solar, about 
1/6 Z Q jKobulnickv et al Jl999t) . The starburst nature of this galaxy finds a direct confirmation in the supernovae observed. Two 
SNe have exploded in less than a cent ury, in 1 895 a nd in 1972, which is a fairly high number given that NGC 5253 is a small 
galaxy: its estimated mass is 6 10 9 M (Welch 1970, scaled to the adopted distance). 

As NGC 5253 has been observed by many authors oyer a wide spectral range for many decades, a precise picture of the galaxy 
can be built today (see e.g. lCaldwell & Phillipslfl989l iMartin & Kennicuttlll995UCalzetti et al.ll 19971 1 1999th We will highlight 
her e only those observat ions that are relevant to our present work. 

iRieke & Low! (1 19721) observed NGC 5253 at 10 and 20 yum as part of their seminal work on infrared photometry of galaxies. 

iGlassI ill 9731) ' and Moorwood & Glassl (1982) also obtained n ear and m id-infrared observations and attempted to constrain the 

spectral energy distribution (SED) in the optical and infrared. iRoche et alJ dl 99 lh . based on their mid infrared spectrum which 

shows strong [SIV] but undetected [Nell], classified NGC 5253, together with II Zw 40, as a high excitation HII galaxy char- 

acterized by a young starburst and by the presence of very massive stars. A more complete mid-IR spectrum was published by 

ICrowther et alJ d 19991) with ISOSWS, and shows that the spectrum is d ominated by fSIVl, [ Nellll and [SHI], a common feature 

of high excitation low-metallicit y HII regions obser ved with ISO (e.g. Contu rsi et alJl20 00 in the SMC). Finally, NGC 5253 is 

detected in all four IRAS bands JMoshiret al.ll990l) . 

1 J 71 

van den Bergh ( 1980) and interpreted as the remnant of a powerful episode 
of star formation. M ore recently high spatial resolution observations have allowed to focus the attention on a few peculiarities. 
ICalzetti et alJ i 19971 hereafter C97) studied the properties of the brightest optical clusters detected in the HST images. They use 
broad and narrow band colors and find that all clusters are young with ages from few to a few tens of Myr. The youngest cluster 
is obscured by about 9 mag in V and in general there is anti-correlation between the presence of dust, that is observed all over 
the central part of the galaxy, and t he age of the clu sters, i.e. the younger the clusters, the high er the extinct i on. Th e centimetric 
radio spectrum is entirely thermal dBeck et all 19961) . which is highly unusual for a galaxy, and Turner al al. (2000) detect a very 
compact radio source, 1 -2 pc in size, which is very bright at 1 .3-2 cm. They attribute this source to a nebula ionized by about 4000 

07 stars. The ionizing flux derived from the radio is about two orders of magnitude larger than indicated by the Ha emission. 

■ 1 j 1 ° 

Goriian et al. (2001) report their Keck telescope detection of a bright infrared source at 11.7 and 18.7 /mi, which they identify 
with the compact radio nebula. 

The paper is organized as follows: in Sect. |2] we describe the new observations on which this work is based. In Sect. [3] we 
discuss the optical colors of the individual bright clusters in the galaxy and exemplify the impact of dust on these colors. In Sect. 
|4]we build and model the spectral energy distribution of the brightest source in the galaxy and show that it can be represented by 
a dust-enshrouded SSC model. The conclusions of our work are summarized in Sect. [5] 

2. Observations 

We have observed NGC 5253 with ISAAC at the ESO-VLT in the Ks and L' near-infrared (NIR) broad bands on April 19, 

2002. ISAAC is equipped with a Hawaii 1024 x 1024 Rockwell detector for the short wavelengths (1-2.5 /mi) and with an 

Aladdin SBRC 1024 x 1024 detector for the long wavelengths (2.5-5 /mi), the scale is 0.148 and 0.071 "/pix respectively. Both 

observations were obtained under seeing better than 0.4". The total on source integration time was 5 minutes in Ks and about 

30 minutes in L' . We have used the usual nodding technique to subtract the background in the Ks band and the nod and chop 

in L' . The images have been reduced both with the ESO package Eclipse and with IRAF, the two independent techniques gave 

perfectly consistent results. The photometric calibration relied on the observation of one photometric standard star. Since the 

Aladdin detector is affected by a significant non linearity, a correction for this effect has been applied to the L' band observations. 

In both filters we have reached a photometric accuracy better than 3% where most of the error is due to the non-uniformity 

across the field. We have c hec ked the consistency o f our photometry and calibrations using the data available in the literature. 

iMoorwood & Glassl dl982h and lForbes et alJ d 19921) provide in their work aperture photometry of NGC 5253 on apertures of 3, 

6, 9, 12, 36 and 48"diameter in Ks and 7.5, 9 and 12"d iameter in L. Using the same apertures on our calibrated i mages, we find 

values within less than 0.01 mag of Forbes et al. (1992) and 0.02 mag on average of Moorwood & Glass ( 1982) in K. In L, we 

can only compare with lMoorwood & G lass ( 1982), the average difference between the two datasets is less than 0.03 mag even 

though large individual variations are present that can be attributed to centering differences and to the large error bars of the 

I II t 

measures of Moorwood & Glass ( 1982). 

In Fig.^we present our Ks and L' band images, the nomenclature for the clusters detected is as in C97. While all the C97 

clusters are detected in the Ks band, only clusters 1 , 3 and 5 are clearly detected in L' . Furthermore we find in L' a relatively bright 

source at about 2.4 arcsec NE of cluster 5 which is not among the bright clusters studied in the optical by C97 and other authors. 

Looking in the HST images we found a faint optical counterpart for this source. The same source is detected in our Ks image. 
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Fig. 1. NGC 5253 in Ks (left) and L' (right) bands. North is up and East is left. Coordinates refer to the J2000 equinox. The field 
of view of the image is 40". The intensity scale on the Ks band has been chosen to clearly show the brightest clusters. The actual 
image is much deeper than this display. The clusters are labeled as in C97. As the numbers in C97 correspond to a brightness 
scale, these images clearly show the dramatic changes occurring in the spectral energy distribution of most clusters in NGC 5253. 
Cluster 5a is the relatively bright object in the L' image North-East of Cluster 5. Note that the intensity levels have been adjusted 
in the L' image to show clusters 1 and 3, but that cluster 5 is a point source in L' as well. Wavelet filtering of the L band image 
reveals that cluster 4 is responsible for the extension seen SW of cluster 5, though the statistical significance of this detection 
is too low to allow a reliable photometric measurement. At an even lower significance level, the same technique indicates the 
possible dection of a source around a = 13 h 39 m 55.6 s and 5 = -31°38'39", but at this very low level, this could very well be a 
false detection. 

Table 1. Aperture photometry obtained with a radius of 0.7"on the main clusters detected in NGC 5253. Typical errors are less 
than 0.03 mag. 





K s 


L' 


1 


14.23 


13.40 


2 


14.60 




3 


14.51 


13.67 


4 


16.10 




5 


13.46 


9.07 


5a 


15.54 


12.51 


6 


15.64 





We will refer this "new" source as 5a. Very faint extended emission, possibly resolved in 2 sources, is observable in L' to the SW 
of cluster 5 in the area of cluster 4 that however cannot be clearly identified. In Table^we present the aperture photometry of the 
clusters studied by C97. For the ISAAC photometry, we did not apply aperture corrections as with the NIR PSFs, the apertures 
used include > 99% of the flux. In addition to those 6 clusters we report the magnitudes of cluster 5a. The photometry for all 
clusters has been obtained with the same method used by C97 for consistency but with a slightly larger aperture to compensate 
for the larger Point Spread Function (PSF) of our NIR observations. In the Ks band the galaxy background has been measured 
on an annulus with radius 0.8"around each source as done by C97 and a width of 0.7 ". In L' no other source than those listed in 
Table^has been firmly detected, neither is any diffuse emission in the galaxy. Cluster 2 and 3 are separated by less than 1 "so 
that their Ks magnitudes are contaminated by each other. 

NGC 5253 has been observed at the SEST with SIMBA on the nights of 7 and 8 of June 2002. SIMBA is a 37 channels 
bolometer array sensitive in the 1.2 mm band. The beam size at the SEST is about 11". The background subtraction is achieved 
by fast scanning of the source through the field of view both in azimuth and altitude. Our maps had a size of 480" x 240". The 
atmosphere opacity r was very good during all observations with values of about 0. 1 . The data have been reduced using the 
package MOPSI developed by Robert Zylka. As a flux reference we have observed the planet Uranus. The total telescope time 
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Table 2. ISOCAM photometry of NGC 5253 



Filter 


FWHM (") 


A mt (jim) 


AA (//m) 


Flux density (mJy) 


LW1 


2.3 


4.5 


1.0 


125±30 


LW10 


3.3 


12.0 


7.0 


1960±300 


LW9 


5.1 


14.9 


2.0 


3900±400 



on NGC 5253 has been about 3 hours which gave a lcr detection limit of 4 mJy. The SIMBA image shows a single unresolved 
source whose flux is 114 + 4 mJy. 

We complement these observations with the analysis of archival data from the mid-IR camera of the Infrared Space 
Observatory, ISOCAM 1 . The interest of ISOCAM data is two-fold: (1) it has a much higher sensitivity than ground-based 
instruments operating in the same range, and (2) with a spatial resolution of a few arcseconds, it helps us build a more pre- 
cise spectral energy distribution for the main cluster of NGC 5253 (see Sect.@}. The data analyzed here were part of program 
LMETC ALE H ARO _A and comprise the identification numbers 62500784, 62500881, and 62500980. We analyzed the observa- 
tions taken with filters LWl-[4.5//m] to serve as a check on our L' observation; LW9-[15yum] to understand the discrepancies 
found when compiling 20/itn photometry (Sect.0}, and LW10-[12yL/m] as this filter is equivalent to the IRAS 12/im filter and 
thus tells us how much of the IRAS flux can be attributed to the central source. All ISOCAM observations were performed 
with the smallest pixel field-of-view, i.e. l'/5, in a "beam-switching" mode, i.e. a neighboring empty field is first observed, then 
the telescope slews to the source. The data were analyzed with CIA V5.0 following the standard procedur es. In all 3 bands the 
source is point-like, which means that the transient-correction scheme used is not adapted (see Coulais & Abergel 2000 for a 
discussion). As a result, our fluxes could be lower limits to the actual flux, though at most by 30-50%, given that a very long 
integration is always performed on the source. The errors quoted in Tableware compound errors, i.e. including photometric and 
transient errors, they are not statistical lcr. Fluxes in Table|2]were obtained with aperture photometry correcting for the missed 
fraction of the PSF, as the ISOCAM PSF has more pronounced low-level wings than the ISAAC one. 

3. Colors of the Bright Clusters 

We have combined our Ks and L' observations with those of C97 to build color-color diagrams of the optically brightest clusters 
in NGC 5253. For source 5a we measured the optical magnitudes on the HST images using the same technique as in C97. To 
ensure consistency between our measures we extracted the optical photometry for all clusters from the HST images, the average 
difference between our values and those of C97 is about 0.02 mag. The F547M-F814W HST color has been converted to V-I 
using the standard photometric method and a linear fit on the data provided in the HST WFPC2 photometry cookbook 2 . It must 
be noted that while the F547M magnitudes are quite close to V, there is a difference of more than 1 mag. between F814W and 
I and a large color term. In Fig.|2]we plot the obs ervations (solid circl es) in V-Ks/V-I and Ks-L'/V-I color-color diagrams. As a 
reference we also plot the output of Starburst99 (iLeitherer et al.lfl999ll - black solid line - calculated for an instantaneous burst 
with a metallicity 1/5 solar and standard IMF (2.35, 1-100 M ). All observations departs significantly from the model showing a 
color excess in all bands, the photometric errors on the colors being 0.1 mag at most for the faintest sources. A warning note is 
due at this point of the discussion. It is well known that sub -solar single stellar population models fail in reproducing the correct 
blue to red supergiants star ratio jLanger & Maederll 19951) and that this has important consequences for the predicted colors 
JOrigliaetalJll 999). For this reason, interpretation of Fig.|2]is not free from ambiguities. The problem however only affects the 
objects in the supergiants phase, cluster 1 and 6 according to the ages derived by C97. In these cases, the red excess could be 
artificial and partially or fully due to the lack of red supergiants in the model. 

To understand the effect that dust can have on the clusters' colors, we examine the color corrections induced by dust emission 
and extinction (a third, more complex effect, i.e. scattering, will be discussed later on). From the colors predicted by the model 
for an age of 1 Myr, the square symbols connected by a continuous line show the effect of adding to the stellar emission an 
increasing fraction of emission from dust at 1000 K (treated as a modified blackbody with an emissivity law in l/A). The square 
symbols indicate the fraction of Ks flux due to dust, with the dust contribution doubling at each step to represent 0, 50, 67, 80, 
88, 94% etc. of the flux. From cluster 5 we plot a screen extinction vector corresponding to a correction for Ay = 2, an open 
triagle indicates the value corrected for Ay = 1. We have used high resolution (7? = 8 10 4 ) optical spectra obtained with UVES 
at the VLT and still unpublished, to measure the extinction from the Balmer decrement. Using a screen extinction model, we 
have obtained Ay = 2.65 and Ay = 2.95 for cluster 1 and 5 respectively. Dust may not present itself only as a screen but could 
also be mixed with the stars, although this geometry may not be too plausible in the case of clusters. In the mixed case, the 
attenuation factor is (l-e~ Tl )/Ti. The color correction terms that result saturate at high optical depth to values dictated by the 
extinction curve that are too small to account for the observed clusters' colors. The screen extinction vector thus represents the 
maximum extinction effect. 



1 The ISO archive is accessible at http: //www. iso . vilspa. esa. es/ 

2 http : /www. stsci .edu/instruments/w£pc2/Wfpcc2_phot/wfpc2_cookbook. html/ 
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Both color-color diagrams give consistent results. The simple examination of Fig. |2 shows that the infrared excess can be 
explained by the combined effect of extinction and emission by warm dust. There is a degeneracy between age, extinction and 
dust fraction which prevents from accurately defining the properties of the clusters. There are however a few important facts: (1) 
an extinction of at least 1-2 mag in V affects all clusters with the only exception of cluster 4; (2) a dust contribution in Ks above 
30% is required to bring the colors of all clusters close to the model prediction (this number is obtained by "sliding" the cluster 
points along the extinction vector till they intersect the dust emission line); (3) clusters 5 and 5a display the most extreme color 
excess: they require more than 2 mag of extinction and a dust contribution in Ks above 70%. The properties of cluster 5 will be 
the main subject of the analysis presented in the following. 

It must be noted that the broad band magnitudes and colors of young star forming regions are usually strongly contaminated 
by the presence of emission lines (the gas continuum emission is taken into account by the Starburst99 model). This e ffect is 
more evident in V where the contribution of the emission lines can easily exceed 0.5 mag jThuan et alJl997|lVanzil2003l) . In our 
case, however, the HST filters F547M and F8 14W are free from bright emission lines so that the V-I color derived from them can 
be considered as representative of the continuum. The only band that can be significantly affected by the presence of emission 
lines is, in our case, Ks so that the V-Ks colors must be considered as upper limits. Since the main emission line present in the 
Ks band is Bry, the effect can be quantified using the Ha fluxes measured by C97 and a standard value for the Ha/Bry ratio. For 
cluster 5, even assuming an extinction A v = 9 as derived by C97, we obtain a correction of only 0.08 mag. The other clusters 
have equivalent widths of Ha that are considerably smaller so that the correction for those cases can be considered negligible. 
Another effect that could potentially affect the V-Ks color is the different angular resolution of the HST and VLT images which 
could be only partially compensated by the different apertures used to extract the photometry. 

Anticipating the results presented in Sect.|4l we draw the reader's attention on the apparent disagreement between the rather 
small visual extinction that can be derived from Fig.[2]and that derived by C97. The reason is that the extinction vector drawn on 
Fig-El as is generally the ca se in such g raphs, only takes into account the extinction effect, and not the scattering by dust grains. 
It has long been known (e.g. lWitt et alJl992i) that while extinction reddens the intrinsic spectrum, scattering makes it bluer. As a 
net result, when scattering is taken into account, the magnitude of the extinction vector decreases for the same Ay. As C97 took 
scattering into account in their estimation of extinction, they obtain a larger value than in the simple extinction-only case plotted 
in Fig|2l On the same figure, we show the color evolution of the model derived in Sect.|4l with an increasing Ty (open circles, 
each circle marks an increase by 1 in the V optical depth, up to ry = 7). This model completely incorporates the extinction, 
scattering and emission effects of dust. Because it incorporates both extinction and emission from dust, the orientation of the 
color correction vector is a combination of those of the emission vector (squares) and screen extinction (triangle) vector. But 
since it also accounts for scattering, the amplitude of the color correction for one mag. of optical depth is smaller than in the pure 
screen extinction case. As will be seen in Sect.|4l our model only partially reproduces the L' band flux, hence the disagreement 
in the Ks-L' color. 

All clusters are barely resolved in our images. In particular in Ks where all field stars display a uniform PSF with FWHM of 
0740, clusters 1 and 6 have a FWHM of 0750, cluster 5 of 0746, cluster 3 of 0752, and cluster 2 of 0760, though this last value 
is made uncertain by the proximity of cluster 3. Cluster 4 is too faint to measure a meaningfull FWHM. Deconvolving with the 
stellar PSF we obtain sizes in the range 0722 - 0745 (3.5-7.2pc) that are fully consistent with those measured by C97. In the L' 
image there are no stars visible so that we have to rely on the standard stars which have a FWHM of 074. Cluster 5, the only one 
for which we can reliably measure the FWHM in L', has a FWHM of 0747 consistent with the value measured in Ks. This means 
that the dust emission is very compact in Ks and L' and that it is confined to the optically emitting region. 

4. Modeling the SED 

As it is evident from figure[2] cluster5 stands out amo ng the popu lation of bright stellar clusters in NGC 5253. Since cluster5 sits 
at the position of the thermal radio nebula of iTurner al all dioOOt) . we can assume that it is indeed the optical counterpart of the 
radio nebula. The fact that optical studies show it to suffer from strong extinction strengthens this hypothesis. Thus in this section 
we first compile observations to build the full spectral energy distribution of cluster 5, from the op tical to the radio, and then we 
attempt to model this SED using the same approach as in SBS 033 5-052 JPlante & Sauvaeel2 002). i.e. modeling the stellar light 
transfer in the dust envelope with DUSTY Jlvezic&Elitzuill997l) . 

4.1. The SED of cluster 5 

From the optical to the L' band our task is simple since the spatial resolution of HST (C97) or VLT (this work) allows to separate 
the cluster from its neighbors and to measure accurately its fluxes. Our optical photometry of cluster 5 agrees within 0.06 mag 
with C97, the difference is mostly due to the background subtraction. 

In the mid-IR range, fl uxes come from a variet y of works, from the early ones with photometers of iRieke & Low! Jl972l). 
i Moorwood & Glass! i ll 982b . and lFrogel etaD (ir982). to state-of-the-art cameras on large ground-based telescopes jGoriian et alJ 
l200ll) and space observatories, namely IRAS and ISO. At the shortest infrared wavelengths, the ISOCAM 4.5 /zm flux of 
125+30mJy compares well with the 77 mJy measured in the L'(3.78 /urn) image, once taken into account the strong rise from Ks 
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0.5 1 1.5 

V-I 

Fig. 2. V-K/V-I and K-L/V-I color-color diagrams of the optically brightest clusters in NGC 5253. The solid black lines represent 
the output of Starburst99 for a instantaneous burst of Z = 1/5Z . The solid lines with square points denote the change on the 
lMyr old point introduced by an increasing amount of dust at 1000K. The arrows indicate the effect of 2 mag of visual extinction 
on the colors of cluster 5. The dashed line with open circles shows how a more proper treatment of all dust effects, extinction, 
scattering and emission, can affect the optical colors. Starting from the 1 Myr dust free Starburst99 model, we compute the colors 
of the same cluster embedded in a cocoon with the dust properties derived in Sect. |4] with the V optical depth increasing by 1 
for each open circle, up to the cluster-5 best-fit value of 7. The photometric errors on the data points are of the same size of the 
symbols or smaller. 

to L' . There is thus no rea son to doubt that all of the ISOCAM 4.5 jum flux comes from cluster 5. Around 10^m the measur ement 
fromlRieke & Low! i 19721) appears rather high compared with the more recent determinations of both Frogel et al. ( 1982) and 
Moorwood & Glass! ill 9821) and we thus choose to exclude it from th e SED. As the SED is seen to increase as the wavelength 
increases, the 2.2 Jy determination at 1 1 .7 fim of lGoriian et al. I fcOOlh can be considered compatible with the rest of the photom- 
etry, although on the high side. The comparison of our LW10 and the IRAS 12/mi flux shows that the ISO flux represents only 
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~75% of the IRAS flux. Given that dust is clearly seen to be widespread in the galaxy, we attribute this discrepancy to diffuse 
dust, to which ISO and ground based telescopes are much less sensitive than IRAS. Diffuse dust, also called cirrus dust in the 
Galactic context, has its emission peak around 100/im. Therefore it must contribute emission to all four IRAS bands, meaning 
that we cannot attribute 100% of the IRAS fluxes at 25, 60 and 100 //m to cluster5. As we do not have high spatial resolution 
measurements at these wavelengths, we assume that the same correction factor applies to them and consider that only 75% of the 
IRAS flux at 25, 60 and 100 /mi originates from cluster 5. 

Moving to the range 15-25 //m, the situation becomes slightly confused. As the 18/mi range is visible from the ground, 
a number of measur ements exist, but they show a rather large scatter, from the small value of 2.9 Jy obtained at 18.7/zm by 
Gorii anet al.l (12001 ) to the large IRAS 25 /mi flux of 9Jy (75% of 12 Jy). Even though the wavelengths are not the same, an 
increase by more than a factor of 4 is quite surprising. In fact we suspect that, although it is backed by some earlier measurements, 
the Gorjian et al. flux underestimates the actual flux in this range. For instance, the point source detected with ISOCAM at 14.9 /mi 
has a flux of 3.9+0.4 Jy, already higher than the 18.7 /mi flux, and very few sources show such a decreasing spectrum from 14.9 to 
18.7 /mi, even when taking into account the second silicate absorption band at 1 8 um. Furthe rmore, and although the photometric 
accuracy on the continuum is probably not very good, the spectrum of Crowth er et al.l (11999) shows that the continuum increases 
slightl y from 15 to 19/imi, with the level at 25 /mi being a good factor o f 2 hig her than at 18/mi. The low Q band flux of 
Riek e&Lowldl972l) at 3.7 Jy could lend some support to the lGoriian et alJ d200ll) measurement however iRieke & Low! dl972l) 
give no detail on the band used for this measurement and thus we cannot use it to constrain the SED. According to Gorjian 
(private communication) the low 18.7 /mi flux could be due to too small a chopping throw, resulting in a n incorrect zero-level, or 
a comp lex airmass correction. We note that a similar measurement in the 18.7 /mi filter was obtained bvlMartin-Hernandez et all 
(2004) giving a higher flux of 5.7 Jy. We will thus discard the 18.7 /mi measurement of Go riian et al.l (120011) from the SED fit. 

Further on in the infrared, we only have IRAS measurements, for which the PSF is almost as large as the galaxy. Again we 
will attribute only 75% of the IRAS flux at 60 and 100 /mi to cluster 5. 

Coming to the submillimeter and millimeter range, we enter a realm where dust continuum emission is not necessarily the 
dominant emission process. Indee d synchrotro n or free-free emission usually dominate longward of a few cm. Here we will 
follow the computations of iTurner et alJ lll997l) . These authors show that longward of 2 cm, and for the compact radio source, 
the flux can reasonably be considered as pure free-free emission. This process has a spectral index of a — -0.1 (in f v oc y"). 
From the 2cm flux they compute the remaining dust continuum emission at 2.6mm, an upper limit of 26+10mJy. This will 
be the longest wavelength included in the SED, and given that it is based on interferometric data we attribut e all this flux to 
cluster5. Using the same method we estimate the free-free contribution in our 1.2 mm observation and in the Lfames et al.ld2002l) 
850 /mi SCUBA observation at respectively 43 and 41 mJy. For those last two measurements however we are left in a situation 
intermedia te between the optical and the IRAS ones: at ~10-20", the PSF could be rather large compared to the object. Judging 
from the llames et al.l d2002l) map, the SCUBA emission comes from the central region of the galaxy, where cluster 5 is located, 
b ut also where most o f the star formation activity, and ionized gas are concentrated (as revealed e.g. by the radio and Ha images 
of lTurneretaiTl 998). This is also the case in our 1.2 mm observations. The association with the ionized gas could simply reflect 
the free-free contribution (~ 20%) but we have no means of ascertaining that. Therefore we will attribute to cluster 5 the full 
free-free-corrected fluxes at 1 .2 mm an d 850/mi, i.e. respectively 71+4 and 151+23 mJy. The errors quoted here are the formal 
errors derived from ours or Ja mes et alJ (120021) photometry. A further source of error comes from the free-free correction, however 
ITurner et alJ dl997l) do not provide enough details for us to estimate the associated error. 

Tableland Fig.[3]summarize all the measurements included in the spectral energy distribution of cluster5. 

Judging from Fig.^one could think that the existence of source 5a brings a further uncertainty to the SED we have just built. 
We do not think that this is the case since already in L' the source is 3.5 mag fainter, and is less red than cluster 5. Furthermore, 
it does not appear to be detected with ground-based telescopes at 10 /mi. 



4.2. Fitting the SED 

Now that the spectral energy distribution has been built, we turn to DUSTY dlvezic & Eli tzurl 19971) to model it. The main feature 
of DUSTY is that it models the radiative transfer of light emitted from a point source inside a spherically symmetric shell of 
dust. Though straightforward to describe, this problem actually involves a rela tively large number o f parameters. The first one 
is obviously the light source. Here, following Sect.[3] we will use Starburst99 (Leither er et alJll9 99) to generate instantaneous 
burst models with ages from 1 to 9Myr at a metallicity of 1/5 Z , close to the one of NGC 5253. Describing the dust shell 
requires many more parameters. The first of these is the dust equilibrium temperature at the inner side of the shell, T m . Due to 
the strategy developed by DUSTY to solve the transfer problem, this temperature sets the physical scale of the model. The radial 
profile of the dust density has to be described. We take it as a broken power law. For N ri d radial zones, this generates {N m <i - 1) 
radius values for the zone limits as the outer radius is fixed in terms of the inner radius, and N m d exponents. To introduce 
some flexibility in the density profile with a limited number of parameters, we chose to run models with two radial zones. In 
principle, DUSTY can include six different dust g rains species. H owever these are mainly variations around three principal types: 
the classical astronomical silicates and graphite (Draine & Lee 1984), and amorphous carbon grains, a possibly more realistic 
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Fig. 3. The spectral energy distribution of cluster 5 in NGC 5253. The 7-30 yum range has been blown up for clarity in the inset. 
Symbols refer to the photometric source, from left to right: black plus signs correspond to the HST photome t ry, gre y diamonds 
to the VLT/I SAAC data, grey triangles to the ISO data, black crosses i n the lO^m range to the lFrogel et alJ (U582) data, black 
circles to the Moor wood & Glassl dl982h data, the black square to the Gorii an et alJ J200ll) 11.7 '/im point, grey plus signs to 
IRAS data, the black cross in the submillimeter range to SCUBA data, the grey square to our SEST/SIMBA observation, and the 
black triangle is the OVRO upper limit. Horizontal error bars describe the filter bandwidths, while vertical error bars show the 
photometric errors, as quoted by the respective authors (see Table[3}. In most cases, these bars are shorter than the symbol size. 
The dashed line is the best fitting DUSTY model. 



representation of interstellar carbon-based dust. With three grains species, the dust composition introduces 2 parameters. Dust 
is also characterized by its size distribution, taken as a power-law, and thus described with three parameters, the minimum and 
maximum grain size, and the exponent of the power law. The standard size distribution, the so-called MRN distribution has 
flmin = 0.005 //m, a mdx = 0.25 /mi (1 jum for graphite) and an exponent q = 3.5 (in the sense n(a) oc g~i. iMafhis et aljfl977h . 
This size distribution is such that, for a 50-50 mix of astronomical silicate and graphite, it reproduces the Galactic interstellar 
extinction curve. Finally, the last parameter describing the shell is its thickness, characterized in DUSTY by the visible optical 
depth ry • 

Therefore to summarize the parameters that describe a DUSTY model in our case are: t s \, the age of the burst, T m the 
temperature on the inner side of the shell, and X gr the abundances of silicates and graphite (the abundance of amorphous 
carbon is X. im = 1 — [Xsi + X„]), a^, a max , and q the parameters of the size distribution, r\ lm the transition radius between the 
two radial zones, fi\ and the exponents of the density power-law (in the sense p oc r~P), and finally r v the optical depth. This 
is a total of 1 1 parameters. As can be seen in Table[3]we have almost 3 times as many photometric data points so the problem is 
clearly over-constrained. However an 1 1 -dimension parameter space is not one that is easy to explore systematically. To identify 
the best-fit model, we choose a strategy where we start from a standard dust shell model and have a restricted set of parameters 
to explore their full space. The best fitting model among this set is then used as a "seed" and a new set of parameters is allowed 
to explore a large part of their respective space. In general, when we start with a new seed model, parameters of the previous 
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Table 3. The compiled spectral energy distribution of Cluster 5. Aperture size are listed when appropriate (i.e. mostly for single- 
beam measurements). A few telescopes and filter bandwidths are missing because we were unable to find the information in the 
original papers. 
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" This measurement is not used as a constraint for the SED fit, see text for details. 
b Only 75% of this flux is used as a constraint for the SED fit, see text for details. 
c Corresponds to the dust-only flux, see text for details. 



exploration are still allowed some margin. This is especially true of t v for purely technical reasons. This is repeated a number of 
times so that all the 1 1 dimensions of the parameter space are explored more than once. 

To compare a model to the observations, the model SED is integrated usi ng the filters' respect i ve tran smis sion curves a nd pho- 
tometric conventions to produce model fluxes. For the photometric data of Moorwood & Glass ( 1982) and Frogel et a 1. (1982), 
we could not find the actual filter transmissions and thus a square bandpass was used. The best fitting model is the one that 
minimizes the quantity: 



JVfflte 



i=l 



f° bS ~f? 



f! 



■obs 



(i) 



The advantage of Eq. is that it gives an equal weight to the small and large flux values, an important fact given the 
very large dynamic range of cluster5's SED (about 6 orders of magnitude), and that it is symmetric, i.e. models that over- or 
under-predict the SED by the same amount have equal weight. 

The best-fitting model is displayed as a dashed line on Figure|3] The search for a better fit in an 1 1 -dimension parameter space 
can be a never-ending one. This model is the best-fitting in the sense that it reproduces the data within the average photometric 
error, with the exception of the F547M filter. In Table|4] we have indicated the range in which individual parameters can vary 
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Table 4. Parameters of the best-fitting dusty model. In the acceptable range, we give an indication of the range of values that, 
while producing a worse fit, still correspond to a model reproducing the data within the photometric uncertainties. 



Parameter 


Value 


Acceptable Range 




1 Myr 


1-2 Myr 
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550-650 K 
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< 0.01 yum 


^max 
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Pi 


0.5 


0-0.5 


TV 
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7-8 



without significantly degrading the fit. Before extracting physical information from the fit, it is worth commenting on some of 
the features discovered as we searched for the best fitting model. 

First, while we have allowed the age to vary between 1 and 9 Myr for a substantial part of the parameter space exploration, 
we have found no intermediate best-fitting model (or "seed" model to use the terminology above) for ages different than 1 or 
2 Myr. This age is younger than that derived fr om the o ptical colors alone, although not by a large amount. It is also smaller than 
the age of the Antennae dust-enshrouded SSC (Gi lbert et al.l2000l) . making it one of the youngest such source. 

Second, the optical depth converges invariably toward a val ue of 7-8. Interestin gly, this is very close to the value measured 
by C97, 9 Vmag. It is also in the same range as that measured bv lGilbert et al.l (2000) for the Antennae SSC. However it does not 
agree with what can be deduced from the color-color diagram presented in Figure|3 the red arrow extending from cluster 5 seems 
to imply only 2 Vmag of extinction. The reason is simple: the effect of dust on the flux in the V, I, Ks and L' bands is much more 
complex than a simple screen extinction on the stellar flux combined with t herma l emission in the infrared bands. Scattering also 
has to be taken into account and, as it is amply demonstrated in IWittet all ill 9921) . it makes the intrinsic spectrum bluer, and thus 
works in the opposite direction of extinction. This very complex effect of dust is exemplified on Fig.|5]with the dashed line: we 
have taken the best fitting model and decreased its visual optical depth (jy) from 7 to almost (i.e. no dust at all). As can be 
seen, the effect on the colors combines that of the simple screen extinction and hot dust contamination, but the color change per 
unit of optical depth is much less than that implied by screen extinction, mostly because of scattering. 

Finally the temperature on the inner side of the shell is also rather stable in all seed models at around 500-600 K, with the 
best-fit value at 570 K. When compared with the results obtained for SBS 0335-052 ( Pla nte & Sauvage1l2002i) . i.e. an optical 
depth of 30 and an internal temperature of 700 K, this shows that the cluster in NGC 5253 is a less extreme case than that of 
SBS 0335-052, though not strikingly different. 

From the actual parameters of the best-fitting model (see Table[4}, we can extract physical parame ters of the embedded 
source and its surrounding cocoon (the equations to perform this can be found in tPlante & Sauvagel2002l) . The best-constrained 
parameter is the luminosity since this is simply the integral of the SED, which is already well described by the observations. The 
luminosity we obtain is 1.2 10 9 L Q . The uncertainty on this value is ~20%. The bolometric luminosity of a nearby galaxy is a very 
difficult quantity to estimate, however one can get a good insight of the importance of cluster 5 by comparing its SED, i.e. Figure|3j 
to the complete SED of the galaxy, readily available on the NASA Extragalactic Database 3 . This shows, quite expectedly, that 
cluster 5 produces nearly all of the IR-Submm luminosity of NGC 5253, but more important, it produces a luminosity equivalent 
to that emerging in the UV-visible region. Therefore this shows that a single source, of very modest physical size, can generate 
as much energy as the complete galaxy. In other terms, it means that even when we consider the global SED of the galaxy, the 
origin of the energy can be a very strong function of the wavelength, i.e. from the widespread distribution of stars over the whole 
galaxy body in the UV-visible, to a single compact source in the IR-Submm. Depending on the representativity of NGC 5253 
with respect to other galaxies, this fact co uld ha ve a significant impac t on all attempts to model the evolution of global SED 
(e.g. ICharv & Elbazll200ll bale & Heloull2002t IC hapman et alj E)03). Finally, to compare clusters to one anothe r, it is often 
conve nient to express the luminosity as a number of equivalent 07V stars (see the properties of this "standard" in lVacca et alJ 
119961) . We find that the bolometric luminosity of cluste r 5 requir es 4700 of these stars, well in the range of SSC's stellar content 
and quite close to the radio estimate of 4000 ( Tu rner al all2 000). 

Using this luminosity and Starburst99 we can obtain the mass of stars present in the cluster. This is 8.2 10 5 M , very close 
to the estimate of C97 (1O 6 M ). This mass obviously depends on the IMF used in the model, in this case a Salpeter IMF with 
masses between 1 and 100 M , and would increase to 2.1 1O 6 M if the mass distribution extends down to 0.1 M . However, 
assuming a single slope value over the whole 0. 1-100 M range is probably not justified (see for instance the detailed study of 



3 see http://nedwww.ipac.caltech.edu/ 
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iRroupa et alll993"l) and usually overestimates the stellar mass. Using a slope 1 .25 between 0. 1 and 1 M G JScalol 19861) . the mass 
of the cluster becomes 1.2 10 6 M . This mass is well in the range observed for SSCs, although compared to SBS 0335-052 or the 
Antennae, it appears on the low side (see Plante & Sauvage 2002; Gilbert et al. 2000). 

Concerning the geometry of the cocoon around the cluster, we find that it extends out to a radius of ~140pc, or approximately 
9". This value strongly depends on T m and it is smaller for larger values of the inner temperature. We point out that this is the 
predicted size of the dust cloud inside which the cluster should reside, but that depending on the observing wavelength, the 
angular size of the source will vary, since the dust temperature strongly decreases with radius. Even though, this is rather large 
compared to the size of the central star formation region of the galaxy (see Fig.Q. This can be resolved in a number of ways. 
(1) Cluster 5 is only located in projection in the central region, but it is in fact behind, or less likely in front of the star forming 
region, or (2) cluster 5 is indeed part of the star forming region, and the dust cloud we model here encompasses it all. In that 
second case, which is more likely, all the central clusters should suffer from some extinction, which is the case, and all these 
provide some heating to the dust, so that our SED for cluster 5 is overestimated. Correcting this is impossible, but would lead 
to a lower luminosity for the central cluster and a smaller stellar mass. We note though, that the high resolution observations in 
the 10-20/mi region show that only cluster 5 contributes infrared emission. This means that the heating provided by the other 
clusters is likely small, hence the correction on the model result should be small as well. It is interesting to note that the predicted 
size should be close to the observed siz e in the submillimeter range, since in that range we are most sensitive to cold dust. Indeed 
the SCUBA map of ljames et al 1 ll2002l) shows a source with a compatible angular size (25-30" diameter at the 6cr level). 

The most striking properties of the dust cocoon, as found from the best-fitting model, are directly connected to the dust 
grains. We find that the exponent of the grain size distribution is 2.5, shallower than in the MRN distribution, i.e. it biases the 
distribution toward large r grains compared to the standard interstellar medium, as was found for SBS 0335-052. As noted in 
IPlante & Sauvage! J2002I) this app ears to be common in high density regions exposed to strong radiation fields (see for instance 
the work of iMaiolino et alJl200ll) . With sizes between 0.005 /im and 0.5 /mi, we do not find in NGC 5253 the lack of small 
grains that was observed in SBS 0335-052 where sizes ranged between 0.2 and 1 fim. Finally the chemical composition of the 
dust is very different from that found in the general ISM of our galaxy: we find that for the best-fit model, the dust composition 
requires no silicate component, being made of 84% graphite and 16% amorphous carbon. This is an extreme compositio n, even 
more severe than in the SBS 0335-052 cases where the b est fit model required 23% silicates dPlante & Sauvagell2002l) . It first 
finds some observational support in the SWS spectrum of Crowther et al 1 lll999l) : even though independent measurements show 
significant extinction toward cluster 5, the spectrum does not show any absorption at the location of the strongest silicate band, 
9.7 fim. We have then explored the deterioration of the fit if we intentionnaly add silicates to the dust mixture, while keeping 
the rest of the model constant. We observe that the figure of merit of the fit degrades gradually, being twice as large as the best 
fit value for a 25% abundance of silicates. However the worst problem is that since our model requires hot grains toward the 
center of the structure to fit the NIR-MIR emission, as soon as the fraction of silicates reaches 10%, we start to observe emission 
features in the model, which are absent from the SWS spectrum for instance. We therefore feel that indeed the dust in the cocoon 
presents a severe lack of silicates. 

We find that the density profile along the shell radius is shallower than in SBS 0335-052. Here the shell has a flat profile 
over the first 2-5% of its radius and then falls with an exponent of 1/2. This profile, combined with the dust composition and 
size distribution, leads to a total dust mass in the envelope of 1.5 10 5 M o . The dust mass is directly proportional to the optical 
depth. It obviously also depends on the chemical composition, dust size distribution and density profile, all three parameters also 
infl uencing the optical depth. We find that the acceptable range of values can lead to a variation of the dust mass by a factor of 
~2. lMeier et al. I (120021) estimate that ~ 10 7 M of H2 is present in the central region of the galaxy. Not all this gas is associated 
to the star formation episode: if we follow the estimation of the same authors, only 5 10 5 M of H2 are causally connected to the 
current starburst, which obviously includes cluster 5. This would imply a molecular gas to dust mass ratio of ~3-4. This is rather 
high, even considering the lack good statistics on this particular ratio. It is however not possible to elaborate on this particular 
issue, given the facts that (1) the molecular mass comes from an interferometric measurement, (2) NGC 5253 is a low-metallicity 
star-forming galaxy where the CO-to-H2 conversion ratio can be quite different from the standard value and (3) the determination 
of which fraction of the H2 gas mass is actually linked with the current burst is a difficult matt er. We also note here t hat the 
mass of dust would be decreased if DUSTY allowed for temperature fluctuations in the dust (see lPlante & Sauvagel20 02). If we 
consider that the dust mass in the cocoon is a good measure of the total dust mass in the galaxy (the cocoon provides almost all 
the IR-Submm luminosity), then the molecular gas to dust mass ratio in the central region of the galaxy becomes ~70, i.e. a more 
reasonable value. 

Finally, we note that in a recent investigation of the ionized nebula associated with cluster 5. lTurner et al.N2003l) consider that 
none of the optical clusters, those detected in the HST data, coincide with the near-infrared source. This rests on a displacement 
between cluster 5 in the HST data and the NIR main source of 073 + 071, which we also confirm with our data. Does this mean 
that indeed the two objects, cluster 5, and the infrared source are different? We think this is too small an offset to guarantee that 
indeed we have two independent sources here. The visual location of the source could for instan ce be affected by n on-uniform 
extinction on the face of the cluster. To be complete however we investigated the consequences iTurner et alJ J2003I) statement, 
which basically result in the HST data dropping out of the fit. The main consequence is that the optical depth has to increase to 
values in the range 15-20, which in turn increases the dust mass required in the enveloppe. We therefore conclude that the model 
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we have presented here represents at least a lower limit (in terms of the mass involved) to the actual situation in NGC 5253. To be 
able to state whether we are dealing with a single source from the visible to the NIR or not would require high-resolution images 
covering the V to K band, in order to "follow" the source's morphology and position with wavelength. 

5. Conclusions 

We can summarize the conclusions of our work as follows: 

- We have obtained very high quality NIR data of NGC 5253 in the Ks and L' bands. The optically bright clusters detected 
all show a red excess that can be explained with the combined effect of extinction and dust emission. This however must 
be considered cautiously for the clusters in the red supergiant phase due to the uncertainty introduced on the models by the 
metallicity. 

- One single cluster dominates the galaxy output beyond 3 //m. We have detected NGC 5253 at 1.2 mm, we attribute the flux 
observed at this wavelength to the IR bright cluster known as cluster 5. 

- Using our observations, archive data and data available in the literature we have built a complete spectral energy distribution 
from the optical to the millimeter, doing this we find further evidence that one single cluster is responsible for most of the IR 
flux which is comparable to the optical flux of the whole galaxy. 

We have then used the model DUSTY to fit the SED and derive the parameters of the IR cluster. We have found that: 

- The IR cluster has a stellar mass of about 0.82-2.1 10 6 M Q depending on the IMF used, and a luminosity equivalent to 4700 
07V stars. This cluster is obscured by 7 mag of optical extinction. 

- The dust in which the cluster is embedded shows the now familiar bias toward shallower size distribution than the MRN one, 
and strikingly lacks silicates, a composition that can explain the absence of a significant silicate absorption band. Adding 
silicates to the best fit model, to obtain a dust composition that resembles the classical one, lead to emission features in the 
model spectrum as the NIR part of the SED requires relatively hot grains. None of these features are observed. 

- The total dust mass necessary to reproduce the optical to submm SED of cluster5 is 1.5 10 5 M o . This is a relatively small 
amount intrinsically, although when compared to the molecular mass observed to be involved in the current episode of star 
formation, it becomes a substantial fraction of the ISM. 

Finally we wish to stress again that NGC 5253 provides a striking example of an object where one part of the spectral energy 
distribution, i.e. the IR-submm one, is completely decoupled from the other, the UV-optical one. The existence of such objects, 
if not a space oddity, would cast some doubt on our ability to predict the global spectral energy distributions of galaxies as they 
evolve in time. 
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